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ABSTRACT 

Context. Linear rotor molecules such as CO, HCO^ and HCN are important probes of star-forming gas. For these species, tempera- 
tures of <50 K are sufficient to produce emission lines that are observable from the ground at (sub)millimeter wavelengths. Molecular 
gas in the environment of massive protostellar objects, however, is known to reach temperatures of several hundred K. To probe this, 
space-based far-infrared observations are required. 

Aims. We aim to reveal the gas energetics in the circumstellar environment of the prototypical high-mass protostellar object 
AFGL2591. 

Methods. Rotational spectral line signatures of CO species, HCO^, CS, HCN and HNC from a 490-1240 GHz survey with Her- 
schelfHWl, complemented by ground-based JCMT and IRAM 30m spectra, cover transitions in the energy range (E„^lk) between 
5K and ~300K. Selected frequency settings in the highest frequency HIFI bands (up to 1850 GHz) extend this range to 750 K for 
'^C'^0. The resolved spectral line profiles are used to separate and study various kinematic components. Observed line intensities are 
compared with a numerical model that calculates excitation balance and radiative transfer based on spherical geometry. 
Results. The line profiles show two emission components, the widest and bluest of which is attributed to an approaching outflow 
and the other to the envelope. We find evidence for progressively more redshifted and wider line profiles from the envelope gas 
with increasing energy level. This trend is qualitatively explained by residual outflow contribution picked up in the systematically 
decreasing beam size. Integrated line intensities for each species decrease as Eap/k increases from <50 to ~700 K. The H2 density and 
temperature of the outflow gas are constrained to ~ 10^-10* cm"' and 60-200 K. In addition, we derive a temperature between 9 and 
17 K and A'(H2) ~ 3x10"' cm"^ for a known foreground cloud seen in absorption, and A'(H2) < 10" cm"- for a second foreground 
component. 

Conclusions. Our spherical envelope model systematically underproduces observed line emission at E^j,/k > 150K for all species. 
This indicates that warm gas should be added to the model and that the model's geometry should provide low optical depth pathways 
for line emission from this warm gas to escape, for example in the form of UV heated outflow cavity walls viewed at a favorable 
inclination angle. Physical and chemical conditions derived for the outflow gas are similar to those in the protostellar envelope, 
possibly indicating that the modest velocity (<IOkms"') outflow component consists of recently swept-up gas. 

Key words, stars: formation - astrochemistry - ISM: molecules - ISM: individual objects: AFGL2591 - ISM: clouds 



1. Introduction 



: Because young high-mass stars (> 8 Mq) reach the hydrogen- 
rN burning phase before the accretion phase ends ( jPalla & Stahler| 
^ [1993 ), their formation process differs from less massive proto- 
stellar objects. Several theories exist for the formation of mas- 
sive stars (see review by [Zinnecker & Yorke||2007| l, including 
stellar mergers and so-called 'competitive accretion' (Bonnell 
|et al.pOOT ; Bonnell & Bate 2005 ). However, the accumulating 
observational evidence for coUimated outflows and disk struc- 
tures in massive star-forming regions ([Beuther et al. 2002, Ce- 
2005) ISandell & Wright|[20T0l |Kraus et al.p010(|Wang| 



saroni 



et al.||20l"2 1 supports a picture of anisotropic accretion through 



* Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



a disk structure, analogous to the low-mass star formation pro- 
cess (e.g., Shu et al!]|1987 Larsonl'2003'), at least for luminosi- 
ties up to 10^ Lq- Even with the latter view of high-mass star 
formation, striking differences arise compared to low-mass star 
formation theory: shorter timescales by a factor ~10 and higher 
internal luminosities by several orders of magnitude. Moreover, 
the bulk of the energy produced by massive stars is emitted at 
UV wavelengths, by which the surrounding gaseous material is 
heated, partly ionized, and otherwise chemically influenced. The 
effect of UV radiation on circumstellar gas has been recognized 
for low-mass protostars (e.g., Spaans et al. 1995) 1 and is even 
stronger in the high-mass case ( Stauber et al.'2004). Since a col- 
lapsing protostellar envelope depends on gas-phase molecules 
and atoms to dissipate gravitational energy, knowledge of the 
chemical composition and excitation conditions of atomic and 
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molecular gas is key in the understanding of the high-mass star 
formation process. 

AFGL2591 is a high-mass star-forming region with a cen- 
trally heated molecular envelope. It is located at galactic coor- 
dinates {,b - 78?9, 0?71 in the Cygnus-X region, which has an 
average distance of 1.7 kpc ( |Motte et al.||2007| . While the dis- 
tance to AFGL259 1 itself has been uncertain in the past ( |Wn| 
der Tak et al. 1999 Schneider et al.||2006^ , in this paper we 



adopt a distance of 3.3 kpc as derived from trigonometric par- 
allax measurements of H2O masers by |Rygl et al.| ( |2012| . With 
respect to the distance of 1 kpc often used in previous work (e.g., 
Hasegawa & MitcheTl' 1995 ; Van der Tak et al. 2000) [Van der| 
Tak & Menten 2005 ; Van der Wiel et al. 201 1 ), the new distance 



changes the estimates for the mass of the dominant protostar, 
VLA 3, to ~40Mq and the bolometric luminosity to 2x10^ L© 
(Sanna et al. 2012; Jimenez-Serra et al. 2012). All linear scales 
mentioned in the remainder of this paper are based on a dis- 
tance of 3.3 kpc to AFGL2591, as is the numerical model de- 
scription we use in Sect.|4] AFGL2591 exhibits a powerful bipo- 
lar molecular outflow extending to >r from the central source 
( |Lada et al.|[T984l [Hasegawa & Mitchell||1995| i. At a 3.3 kpc 
distance, this corresponds to a linear scale ~ 1 pc, which is not 
untypical for massive protostellar outflows (Beuther et al. 2002 | l. 
A massive disk has been proposed to exist around source VLA 3, 



but with an angular extent of <1" (Van der Tak et al. 2006 
Wang et aL||201 2 ) the disk will not be resolved in the current 
study, which traces >10" scales. Although interferometric ob- 
servations have revealed multiple sources inside the envelope at 
spacings of <10000 AU ([Van der Tak et al.|1999| [Trinidad et al. 



[2U0T; Benz et al. 2 007[|De Wit et al.|2009t [Sanna et al.|2012i l 

AFGL2591 is isolated at scales of a few parsecs. Such isola- 
tion is exceptional for massive star-forming envelopes, which 
are usually found in clustered environments. The absence of 
complex dynamical and radiative interaction with neighboring 
objects has made AFGL2591 a popular high-mass protostellar 
object to study. 

Most previous modeUng efforts for the large scale envelope 
of AFGL2591 have assumed a spherical morphology with a 
power law density profile ([Van der Tak et al. 1999 2 000^ Doty 
et al. 2002; Stiiuber et al. 2005; Benz et al. 2007; D e Wit et al. 
2009). In these studies, gas temperatures range from ~20K in 
the outskirts of the envelope (~10^ AU from the center) to sev- 
eral hundred K at distances as close as a few hundred AU from 
the protostar. Likewise, Van der Wiel et al.[ ( [201 1 ) employ spher- 
ically symmetric model geometries, but they also turn to flat- 
tened axially symmetric geometries to explain molecular emis- 
sion maps at scales of ~10^-10^AU. Moreover, after sugges- 
tions for the presence of cavities in the envelope carved out by 
outflow motions ( |Van der Tak et al.ll999 1, and observational ev- 
idence in the near-infrared by Preibisc h et al.[ ( [2003] l, detailed ra- 
diative transfer models including outflow cavities are presented 
by [Bruderer et al. ( 2009a|b 2010a) l. They include a chemical 
balance to allow for photodissociation of molecules along the in- 
terface between the envelope and the cavity, and show that a con- 
siderable fraction of the gas is at temperatures exceeding 1000 K. 
This adjusted temperature balance and the possible presence of 
shocked regions obviously affect the excitation conditions of the 
molecular gas. 

In the envelopes of massive star-forming regions, much of 
the gas is well below 50 K, but temperatures also exceed 100 K in 
a significant portion of the gas. In this work, we study this warm 
gas in the large scale envelope of AFGL2591. To study these 
high gas temperatures with ground-based observations, one re- 
lies on low abundance species with relatively complex spectra 



such as CH3OH, SO2 and H2CO. Because H2O is only abundant 
at temperatures above ~100K ( Ceccarelli et al. [[19961 Doty & 
,Neufeld 1997), the most important probes of star-forming gas 
are linear rotor molecules. These species, however, are only ob- 
servable from the ground in transitions that trace either very hot 
(~10^ K for near-infrared H2 lines) or cold gas (<50K for low J 
transitions of asymmetric linear rotors in the (sub)millimeter). 
Space-based observations are therefore necessary to obtain a 
comprehensive picture of the warm gas at a few hundred K. We 
exploit the far-infrared frequency range accessible to the Her- 
schel Space Observatory (Pilbratt et al. 2010) and its Heterodyne 
Instrument for the Far- Infrared (HIFI, De Graauw et al. 2010) to 
observe CO, HCO+, CS, HCN and HNC in high rotational tran- 
sitions with upper level energies ranging from 80 to ~700 K. 

This paper is structured as follows. Section |2] describes 
the observations and processing of spectra obtained with Her- 
ic/ieZ/HIFI from space and of complementary ground-based 
millimeter-wave observations. In Sect. [3] we present line detec- 
tions of CO and isotopologues, HCO+, CS, HCN and HNC, and 
we discuss trends in line position, shape and intensity. The ob- 
served line strengths of the envelope component are compared 
with radiative transfer models in Sect. |4] Section |5] presents a 
derivation of physical conditions in kinematically distinct com- 
ponents. Finally, we discuss the results in Sect.|6]and summarize 
our conclusions in Sect. [7] 



2. Observations 

2.1. 490-1850 GHz Herschel/HIFI spectroscopy 

The Herschel Key Program "Chemical HErschel Surveys of 
Star-forming regions" (CHESS, Ceccarelli et al.|2010| l uses the 
HIFI instrument to perform spectral surveys of various Galac- 
tic star-forming regions. One of the targets is AFGL2591, for 
which an unbiased spectral survey was conducted in the 490- 
1240GHz range (bands la-5a, totaling 18.4 hours of observ- 
ing time), complemented by selected frequency settings between 
1260 and 1850 GHz (bands 5b-7b). Of the total of eight targeted 
frequency settings, this work only uses the five that cover '^CO 
transitions (2.1 hours of observing time). 

The observed frequency bands are listed in Table [T] along 
with observing dates, typical values for the noise level, system 
temperature, and half -power beam width (HPBW). Observations 
were carried out using the dual beam switch mode with off po- 
sitions 3' east and west of the target position: RA - 20''29"2459, 
5 = 40°11'21" (J2000). The specti-al scans (bands la-5a) were 
conducted with the wideband spectrometer (WBS) using a re- 
dundancy factor of 4. Per band, per polarization (H or V), per 
sideband (lower or upper), the scans contain 90-151 datasets, 
each composed of four 1 GHz 'subband' spectra. The pointed 
frequency settings (bands 5b-7b) were also executed in the dual 
beam switch mode, in this case with the fast chop and stabil- 
ity optimization options selected. The spectral resolution of the 
WBS is 1.1 MHz, coiTesponding to 0.66 km s ' at 500 GHz and 
0.18kms-i at 1850 GHz. 

The raw data of the spectral scans were processed with the 
standard pipeline included in HIPE (Ott 2010) version 8.1. Fur- 
ther processing of the individual datasets in the spectral scans 
was done in HIPE, using scripts written by the CHESS team 
( |Kama et al.|2013 1: (1) data inspection and flagging of spectral 
regions with spurious features, (2) subtraction of sinusoid stand- 
ing waves, (3) subtraction of a polynomial baseline, generally of 
order < 3, fitted to line-free regions. Strong spectral lines are 
known to create minor ghost lines in the sideband deconvolution 
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Table 1. Overview of Herschel/KIFI observations of AFGL2591 used in this work. 



Band 


observing date 
(yyyy-mm-dd) 


Herschel obsid 


Freq. range" 
(GHz) 


'obs 

(s) 


rms noise"^ 
(mK) 


d 

sys 

(K) 


HPBW'= 
(") 


77mK^ 

'/ITID 

(%) 


la 


2010-04-11 


1342194483 


483-558 


4592 


30 


74-11 


44-38 


75.5 


lb 


2010-04-12 


1342194528 


555-636 


4643 


29 


8r\l 


38-33 


75.3 


2a 


2010-10-28 


1342207599 


631-722 


9833 


26 


142-11 


34-29 


75.1 


2b 


2010-04-12 


1342194573 


717-800 


6407 


67 


188!^6 


30-27 


74.9 


3a 


2010-10-28 


1342207622 


800-859 


4893 


39 


200^1 


27-25 


74.7 


3b 


2010-04-13 


1342194699 


858-960 


8578 


67 


205+148 


25-22 


74.4 


4a 


2010-05-16 


1342196594 


950-1060 


9137 


157 


356-11'^ 


22-20 


74.1 


4b 


2010-05-11 


1342196428 


1051-1120 


6300 


144 


339!^™ 


20-19 


73.8 


5a 


2010-05-12 


1342196510 


1110-1240 


11931 


147 


qi T+343 


19-17 


63.7 


5b 
6a 
6b 
7a 
7b 


2011-05-11 
2011-04-21 
2010-12-02 

2010- 12-04 

2011- 05-11 


1342220520 
1342219224 
1342210726 
1342210803 
1342220471 


-1267 
-1497 
-1612 
-1727 
-1841 


1380 
1440 
1392 
1575 
1711 


149 
117 
106 
92 
92 


1246^126 
1507^143 

1366±'°4 

1300^43 
1201^5 


17 
14 
13 
12 
12 


63.4 
71.8 
71.1 
70.5 
69.7 



Observations from bands la to 5a are spectral scans, those in bands 5b to 7b are single frequency settings spanning -4 GHz. 
Observing time, including overheads. 
Typical rms noise level in 0.5 MHz channels in T^b units. 
^'^^ Double sideband system temperature. For spectral scans (la-5a): median of the individual datasets, with superscripts and 
subscripts denoting the deviations to the 90* and 10* percentiles, respectively. For pointed frequency settings (5b-7b): average of 
the values for the H- and V-datasets, with superscripts denoting the devia tion to lowest and highe st value. 



Half -power beam width at the relevant frequency, calculated following Roelfsema et al. ( 2012| . The two values for each spectral 
scan indicate the HPBW at the lower and upper edge of the band. 
*0 Typical value of the main beam efficiency in the center of the band. 



process (see below). Therefore, the above steps were repeated 
by flagging not just the spurious features, but also the strong 
spectral lines, resulting in a cleaner secondary spectrum with all 
strong lines masked. The threshold for marking lines as strong 
varied from Ta* =1 K in band 1 to 8 K in band 5a, depending on 
the amount of lines and the noise level in a particular band. 

The spectral scans of the CHESS program were executed 
such that each sky frequency was covered by several settings 
of the local oscillator, which allows for the mathematical dis- 
entanglement of the upper and lower sideband signals that re- 
sult from double sideband heterodyne observations ( jComito &| 
[Schilke 2002). This 'sideband deconvolution' was performed 
with sideband gain ratios fixed to unity. Spectra from the HIFI 
spectral scans (490-1240 GHz) presented in this paper are all 
single sideband, 'unfolded' spectra. 

The correspondence of deconvolved products for the two po- 
larization backends, H and V, is checked and found to agree to 
within 3% in peak intensity for various lines, and line peak posi- 
tions and line shapes generally match to within the noise level 
(cf. Table [T]i. The signal from the two polarizations is aver- 
aged by feeding the datasets from the two backends to the side- 
band deconvolution algorithm simultaneously. One exception is 
band 4b, where we used only the H polarization, since V is con- 
taminated by numerous spurious features and the V-only decon- 
volved spectrum exhibits noise levels of up to four times higher 
than in the H polarization. 

The WBS spectra from the pointed frequency settings of CO 
transitions above 1240 GHz (bands 5b, 6a, 6b, 7a and 7b) are 
processed with HIPE version 7. Subsequent removal of poly- 
nomial baselines of order <4 and sine-shaped standing waves is 



applied in the CLASS packag^ With a single frequency setting, 
sideband deconvolution is impossible, so the high-frequency CO 
lines are double sideband spectra. Like in most spectral scans 
described above, the average of the two polarizations is used, 
with the exception of CO 11-10 in band 5b. For this line, the 
spectral line shape shows discrepancies between H and V. Since 
the baseline rms noise level is -20% higher in V, we only use 
the H polarization in this case. 

We use the main beam efficiency r]^\„ as determined from 
a beam characterization experiment based on Mars observations 
( Roelfsema et al. 2012 1, to convert the T/ intensity scale to main 
beam intensity: T^\, - {rif-v,dl'nmb)TA' with the forward efficiency 
y/fwd = 0.96 for HIFI. Although ?7,nb is in principle frequency- 
dependent, it varies by < 1 % across an individual band. A more 
significant variation in the value of /7mb is the 14% drop between 
band 4b and 5a (Table[TJ. 

The scaling to main beam intensity is strictly only applica- 
ble to emission regions that are not larger than the main beam of 
the telescope. For the lines studied here, this criterion is largely 
obeyed, so we apply the rj^^, factor to facilitate a direct compar- 
ison to modeled intensities. 

Conservative estimates put the formal uncertainty for the in- 
tensity calibration of HIFI at 9% in bands 1 and 2, at 10% in 
bands 3 and 4, and -14% in bands 5-7 (Roelfsema et al. 2012). 
However, the uncertainty budget is often dominated by stand- 
ing waves and the uncertainty in the sideband ratio, both of 
which are mitigated considerably by our systematic cleaning (all 
bands) and sideband deconvolution of the spectral scans (bands 
la-5a). Based on consistency checks of HIFI line measurements 

' CLASS is part of the GILDAS software developed at IRAM: 
http : //www . iraiti . fr/IRAMFR/GILDAS. 
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over Herschel's operational lifetime, the observed depth of satu- 
rated absorption lines, comparison to ground-based observations 
and cross-comparison with Herschel/SPIRE, realistic values for 
HIFI's absolute calibration uncertainty should be between 5 and 
10%. We adopt a uniform uncertainty of 10% for the line inten- 
sities measured from our HIFI spectra. 

2.2. Ground-based (sub) -millimeter spectroscopy 

Apart from the Herschel spectra, this paper uses ground-based 
observations with the IRAM 30m Telescop^ and the James 
Clerk Maxwell Telescop^(JCMT). 

IRAM 30m observations were performed on July 21, 2010. 
The observations were conducted with the EMIRO and EMIR2 
receivers in parallel (two linear polarizations each) with the 
VESPA autocorrelator backend in the 2.7 mm (~110GHz) and 
1.3 mm (-230 GHz) window^ The C'^^O and '^CO 7=1-0 Hnes 
were observed simultaneously by splitting the backend in two 
subwindows centered on the lines with a frequency sampling of 
78 kHz (-0.2 km s^'). The C"*0 2-1 line was obsei-ved with 
the same frequency sampling (~0.1 kms"'). We used frequency 
switching mode and an integration time of 1 minute on the sky 
(2 minutes after averaging the two polarizations). The weather 
was good with precipitable water vapor of 3 mm, and system 
temperature was 100 K at 2.7 mm and 200 K at 1.3 mm, in Ta* 
scale. A fourth order polynomial baseline has been subtracted. 
The measured noise in the spectra, in T^b units, is 46 mK for the 
2.7 mm spectra and 116mK for the 1.3 mm spectrum after av- 
eraging and folding. We use HPBW of 22" and 11", and main 
beam efficiencies (z/mb) of 78% and 61% for 2.7 mm and 1.3 mm, 
respectively. For the IRAM/EMIR spectra, we take an absolute 
intensity calibration uncertainty of 10%. 

Spectra of rotational transitions of CO species, HCN, HNC, 
HCO+, and CS that fall between 330 and 373 GHz are extra cted 
from JCMT Spectr al Legacy Survey data ( [Plume et al.|2007[|Van| 
|der Wiel et al.|2 01 1 ). The beam size of the JCMT is 14-15" at 
these frequencies, and the spectral resolution used in the sur- 
vey is 1 MHz, corresponding to ~0.8 km s ' . The absolute inten- 
sity calibration for the JCMT measurements is uncertain to 15%. 
A detailed description of the data acquisition and processing is 
given by |Van der Wiel et al.| ( [201 \\ . 



3. Observational results 

3. 1 . Detected emission lines and line shapes 

The spectral scans and targeted frequency settings described in 
Sect. |2] cover high-/ transitions of the linear molecules HCOJ, 
CS, HCN, HNC, CO, '"'CO, C"*0, and C'^O, as listed in TableU 
Other detected spectral lines will be presented in a companion 
paper by Kazmierczak et al.| ( |in preparation! hereafter 'Paper II, ). 

For all but one of the species mentioned above, frequencies 
of rotational transitions are extracted from the Cologne Database 
for Molecular Specti'oscopy (CDMS, Miiller et al. 2005). For 



~ IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
ION (Spain). 

' The James Clerk Maxwell Telescope is operated by the Joint Astron- 
omy Centre on behalf of the Science and Technology Facilities Council 
of the United Kingdom, the Netherlands Organisation for Scientific Re- 
search, and the National Research Council of Canada. 

See |http : //www . iram- institute . org/| for instrument specifica- 
tions. 
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Fig. 1. Spectral line profiles of CO, "CO, C'*0, and C"0. The 1-0 
and 2-1 spectra are from the IRAM 30m telescope, the 3-2 lines are 
from JCMT/HARP-B. All higher-i lines are from Herschel/mFl. For 
clarity, consecutive spectra, each labeled by their upper and lower J- 
level, are offset in the vertical direction. Corresponding rest frequen- 
cies and energy levels are listed in Table |2] The systemic velocity 
(-5.8 km s"') of AFGL2591 is indicated by the dashed vertical line. 
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Fig. 2. As Fig.[T] but for HCO+, CS, HCN, and HNC. All lines are 
measured with Herschel/HIFI, except HCO+, HCN, and HNC 4-3, as 
well as CS 7-6, which are measured with JCMT. 



C'^O we use the JPL frequencies ( [Pickett et aL]|l998|p] We 
search the processed HIFI spectra for line signal near these fre- 
quencies. For weak lines (peak intensity < ScTrms) the spectrum 
is smoothed in steps of a factor 2, until channels become so wide 
that they become comparable to the line width. If the S /N is 
still below 5 at this point, a conservative upper limit to the in- 
tegrated line flux is then set at Scr^^s x (3kms"'), with the last 
factor motivated by the typical line width expected from firmly 
detected lines at lower J. 

Spectral line profiles of transitions of CO, '^CO, C'**0, and 
C'^O are presented in Fig. [l] those of HCO+, CS, HCN, and 
HNC are shown in Fig. |2] (see also Table |2]). If a line needs 
smoothing of the spectrum before its peak intensity is above 
5o"i-ms, Fig.[2]gives the smoothed spectrum. 

All lines with sufficient signal-to-noise show an asymmet- 
ric profile with line wings that are stronger on the blueshifted 
than on the redshifted side. The blueshifted line wing (Visr < 
-lOkms"') is visible in emission in the lowest HIFI transitions 
of CO, and to a lesser degree of HCO^. We attribute the asym- 



^ The JPL database lists an average frequency for each rotational tran- 
sition of C"0, whereas CDMS lists all 15 hyperfine transitions sepa- 





-10 -5 

(kms-i) 

Fig. 3. Examples of Gaussian component fits to HCN 6-5 (top, single 
component) and HCO* 6-5 (bottom, two components). Observed line 
profiles are shown as histograms, fit components are shown as smooth 
solid curves and the sum of the two components in the bottom panel is 
represented by the dashed curve. 



metry on the bl ue si de of the line profile to the approaching out- 
In both CO and HCO+, 



rately, which - at a maximum separation of ~2 MHz 
all C"0 lines observed in this work. 



are blended in 



flow (see Sect. 3.2 1. In both CO and HCO , it becomes less 
pronounced at higher /-levels, but it is nevertheless present up 
to the highest CO transition observed (7=16-15). A separation 
between envelope and outflow emission is also seen in CO 2-1 
spectra by Mitchell et al. ( 1992 ), similar to our CO 3-2 profile. 
In addition, this line wing is known from infrared absorption 
studies of '^CO and '^ CO. In '^CO, the out flow is seen out to 
200k m s'' blueshifted ( [Mitchell et al.|l989|[\^ der Tak et aL 
[T999|. 

Contrary to most bright lines, which peak at the systemic ve- 
locity and have a blueshifted line wing, the peaks of all '^C'^O 
lines (Fig.[T]i are displaced toward redshifted velocities. This is 
attributed to the optical depth and the possible presence of ab- 
sorption components affecting the line shape. Both these effects 
scale with the overall abundance, which is far higher for '^C^^O 
than for any other species studied here. In addition, the 12-1 1 
and 14—13 lines of CS (Fig.[2]i appear asymmetric with excess 
emission on the redshifted side at the ~2-cr level. 

The 11-10 line of HCO^ appears to be split in two separate 
components. The brightest one follows the trend of increasingly 
redder peak velocity with higher energies, while a secondary, 
narrow peak emerges at -8kms"' (5cr detection). The same 
feature, although fainter than for 7=1 1-10, is visible in the 10-9 
and possibly also the 9-8 lines of HCO^ when a single com- 
ponent profile based on lower lines is subtracted. In the other 
species studied here, this additional narrow component is found 
only in C'^O 8-7 at the 2.5-cr level, while it might be hidden 
in the line wing of C'^O 7-6. Based on the firm detection in 
HCO^, we conclude that the feature must originate from dense 
and very warm gas, given that it is brighter at 7=1 1-10 than at 
lower transitions. 
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Fig. 4. Spectral line energy distributions of CO (a), its isotopic variants (b, c, d), HCO^ (e), CS (f), HCN (g) and HNC (h), as a function of 
iiup. Error bars on the observations (in black) represent the quadratic sum of the absolute intensity calibration uncertainty and the formal error on 
the Gaussian fit. Observed intensities represent only the envelope component for those lines fitted with a double Gaussian profile (Table|2|. Since 
the intensity for a line fit by a single Gaussian is likely to be contaminated by outflow contribution as well, gray error bars provide an indication 
of the line intensity when up to half of the total line strength is attributed to the outflow. Downward pointing triangles indicate upper limits for 
undetected lines (see text in Sect. |3.1[ ). Integrated line intensity values from the spherical model are marked by crosses and connected by dashed 
lines. Intensities from ground-based observations (7=1-0, 2-1 and 3-2 for CO species, lowest transitions shown for HCO^, CS, HCN and HNC) 
and the corresponding model points (left of the dotted vertical line) are scaled down by a factor 10, as indicated by the '-=-10' labels. Gray crosses 
accompanied by upward pointing arrows denote cases where the modeled line intensity falls off the vertical scale. 



3.2. Observed spectral line energy distributions 

Integrated line intensities, J TmbdV, are determined by fitting 
Gaussian profiles to the observed spectral lines (Table [2]), using 
the Levenberg-Marquardt fitter in the line analysis module of 
CASsif ^ First, a single Gaussian profile is fitted to each observed 
line. ^ In cases where this yields a residual spectrum with an 
rms noise comparable to or below the local baseline noise level, 
a single Gaussian profile is adopted. If subtraction of the sin- 
gle Gaussian fit leaves a significant residual signal, a fit of two 
blended Gaussian profiles is performed. If this gives a residual 
spectrum with a significantly lower rms than the one resulting 
from the single Gaussian fit, we adopt the double Gaussian fit. 
Examples of cases where a single or a double Gaussian fit is 
appropriate are shown in Fig. [3] Although the CO lines up to 
may accommodate more than two Gaussian components. 



* CASSIS is developed by CESR-UPS/CNRS (http://cassis.cesr.frl. 
' Only in the case of CO 16-15 in band 7b is the Gaussian line fitting 
preceded by subtracting a linear baseline to minimize the influence from 
the nearby OH line in the other sideband. 



the additional components are so faint that they have negligible 
effect on the two primary components. 

Lines that are well represented by a double Gaussian fit are 
CO, '^CO, C'^O and HCO+ up to a certain J level, as well as all 
ground-based line profiles, while a single Gaussian suffices for 
the HIFI lines of C'^O, HCN, HNC, and CS. In general, a lower 
signal-to-noise of a spectral line reduces the chance of finding a 
meaningful secondary Gaussian component. For the lines with 
a double Gaussian profile, we attribute the component with the 
narrowest line shape and reddest Vcentr to the quiescent envelope, 
and the broader blueshifted component to outflow material. For 
lines with a single Gaussian profile, all emission is attributed to 
the envelope. For CO and HCO^, the component separation al- 
lows for an analysis of variations in physical conditions between 



the envelope and the outflow. This is addressed in Sect. 5.2 



Figure|4]shows the J T^^jdV of the envelope components ver- 
sus the energy of the upper level (fiup). For the lines where only 
a single Gaussian fit is used, it is still likely that up to ~50% of 
this emission stems from the outflow component. We use a gray 
errorbar to indicate the intensity level that remains after 0-50% 
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of the single Gaussian intensity is removed from the envelope 
component intensity. The general trend in these spectral line en- 
ergy distributions (SLEDs) is that the line intensities decrease 
with increasing upper level energy. A detailed comparison of 
the observed envelope SLEDs to models is discussed in Sect.|4] 
We emphasize that, although we present various trends as a 
function of E^p, it is important to recognize the additional sys- 
tematic effects of decreasing beam size, and decreasing line op- 
tical depth, simultaneously with increasing E^p for consecutive 
transitions. 



3.3. Observed centroid velocities and line widths 

From the fitted Gaussian profiles, we also obtain measures for 
centroid velocity and line width (Table [2]i. Figure [5] shows the 
trend of fitted values of Vcentr and FWHM for the envelope com- 
ponent, as a function of the upper level energy of the line tran- 
sition. We only plot those species for which at least two lines 
are fitted with a double Gaussian profile. The main form of CO 
is excluded, since its line widths are heavily influenced by op- 
tical depth eff'ects. Results from the HIFI spectral scans of the 
low-mass protostar IRAS 16293, also from the CHESS program, 
show a similar trend as is found in our data of AFGL2591, but 
with Vcentr decreasing with increasing E^p (E. Caux, private com- 
munication). We therefore conclude that the shifting centroid 
velocity over the HIFI bands is not an instrumental eff'ect. 

In the diagram for Vcenti (Fig- |5] top), the two CO iso- 
topologues appear to have a minimum Vcentr (bluest peak) at 
£up~80K (J-5^), with a monotonously increasing trend to- 
ward higher E^^p. The four HCO^ lines also exhibit an over- 
all trend of centroid velocities becoming redder with increasing 
fiup- To explain this systematic shift in centroid velocity, we note 
that increasingly higher J levels do not only probe higher tem- 
peratures, parameterized by E^p, but also progressively smaller 
portions of the envelope, due to the decreasing angular beam size 
(Tables [1] and [2]i. In cases such as the envelope of AFGL2591, 
which is internally heated and enjoys a radially declining tem- 
perature profile, increasing J level and decreasing beam size 
both work in the same direction: the line signal probes pro- 
gressively warmer gas. The observed velocity trend could be 
explained by (i) the very warmest core of the envelope being 
at a slightly different Visi than the surrounding envelope, or (ii) 
colder outer parts of the envelope being (mildly) influenced by 
blueshifted outflow material. The latter explanation is corrobo- 
rated by the '^CO and C"*0 lines up to 3-2, which have redder 
Vcentr than their 5-4 lines observed by HIFI (Fig.[T]). The lines at 
£'up<50K lines are observed at IRAM 30m or JCMT with beam 
sizes of 22" ('^CO and C'^^O 1-0), 11" (C'^^O 2-1) and 15" 
('^CO 3-2 and HCO+ 5-4), afl much smaller than the 38" beam 
applicable for the lowest-/ HIFI lines. 

To test the possible influence of the beam size on the centroid 
velocity value, we take advantage of the I'xl' map of C'**0 2-1 
obtained at the IRAM 30m telescope. When the map is con- 
volved to simulate the largest Herschel beam, the velocity cen- 
troid is measured to move from -5.60 to -5.85 km s"', demon- 
strating that probing a larger portion of the envelope results in 
bluer line profiles. We conclude that the trend of increasing Visr 
with increasing 7-level is likely to be caused by the systemati- 
cally decreasing Herschel beam size (see Fig. [6]). 

In addition, we investigate the trend of FWHM of the lines as 
a function of E^p in Fig.|5](bottom), as measured in the envelope 
component of the Gaussian fits. Lines of '^CO and HCO^ are 
systematically narrower with increasing E^^p and, again, decreas- 
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Fig. 5. Trends of fitted centroid velocity (top) and fitted FWHM 
(bottom) of the envelope component versus £„p. 




Fig. 6. Indication of the half-power beam width of /ferac7je//HIFI, rep- 
resented by black circles: 42" at 500 GHz (left) and 19" at 1100 GHz 
(right). The contour map represents HCO^ 4-3 emission toward 
AFGL2591 measured by the HARP-B instrument at JCMT. The blue- 
and redshifted wings are shown in solid and dashed contours, respec- 
tively, and are defined as in Van der Wiel et al.| ( |201 1) . The lowest con- 
tour and spacing between consecutive contours is 5Kkms"' for both 
components. 



ing beam size. As with the Vcentr trend, the two ground-based 
lines for '^CO form an exception. The systematic shift is less 
clear or even reversed for C'^O. 

To explain systematically decreasing line widths, a turbulent 
core model ( |McKee & Tan|[2003| l could be invoked, in which 
turbulent velocity increases with radial distance due to energy 
dissipating from the warm inner regions into the colder outer en- 

( 2012| l successfully model 



velope. For example, Herpin et al 



an extensive set of H2O lines with such a model description. If 
we adopt a similar model, however, line widths would be ex- 
pected to increase toward lower transitions for all species. Since 
the observed FWHM values for C"*0 between and 150 K and 
'^CO below 50 K show no such relation, we fall back on the 
hypothesis that the outflow material contributes to the 'envelope 
component' line widths. As illustrated in Fig.|6] the smaller tele- 
scope beam at higher frequencies (and higher E^p) picks up less 
emission from the spatially extended outflow. 

In addition, systematic trends of line velocity and line width 
with energy and/or beam size might be explained by an infall 
velocity gradient. This does not seem to be the case, however, 
for AFGL2591, because, even at a spectral resolution of 0.1- 
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1 km none of our line shapes (Figs.[I]|2|i show asymmetries 
that could indicate infall. 

In conclusion, the trends in both Visi and FWHM suggest 
that the 'envelope component' is kinematically affected by the 
outflow material, even after separating the obvious contribution 
of outflow material by means of the double Gaussian fits. This 
implies that our method to separate the outflow from the enve- 
lope is only partially successful. The assumption that the two 
physical components, quiescent envelope and outflow, are fully 
disconnected and do not influence each other dynamically obvi- 
ously constitutes a highly simplified description of reality. 

3.4. Angular size of emission region 

Information about the angular size of an emission region can 
be derived from direct comparison of line intensities measured 
at the same frequency with different telescopes and beam sizes. 
One of the few atmospheric transmission windows that allows 
far-infrared frequencies to be observed from the ground is the 
region around 800 GHz, where this work presents the HCN 9-8 
lines at 797.433 GHz and the CO 7-6 fine at 806.652 GHz in 
a 27" Herschel beam. These two lines were measured toward 
AFGL2591 by |Boonman et al.|pOOT) using the JCMT, with a 
beamsize of 8" (FWHM). The angular size of the emission re- 
gion (0s) is coupled with the measured line intensities through 
its relation with intrinsic line brightness T^, observed main beam 
line brightness r^b, and telescope beam size flbeam, ( [Wilson et al.| 
[2009] l: 



Table 3. Comparison of peak intensities of lines observed by both Her- 
schel (3.5 m) and JCMT (15 m). 



s ^beam 



(1) 



assuming that 9^ is small compared to the telescope beam. Be- 
cause should be independent of the observing facility, 9^ can 
be derived by substituting sets of (r^b, 6'beam) for both JCMT and 
Herschel and equating the results. 

We take the peak T^^ values of CO 7-6, HCN 9-8 and 
HCN (V2=l) 9-8 as measured at the JCMT, with an uncertainty 
of 50% ( Boonman et al.|[2001) l. In Table |3] we compare them 
to the peak intensities of the same lines measured in our Her- 
schellHWl spectrum. For the vibrationally excited HCN, the 
HIFI line strength is an upper limit. Equation ([T} is then used 
to derive limits on the angular size of the emitting region, 0s- In 
the limiting case of point-like emission regions in both beams, 
i.e., 9s^ ^ 0beain^, the line intensity ratio ]CMT/Herschel would 
be 27^/8^ X! 1 1, from simple beam dilution. For larger emission 
regions, the intensity ratio decreases. 

The values for ^s listed in Table [3] imply that the HCN 9-8 
emission originates in a very confined region of the envelope of 
at most 1" in diameter. The CO 7-6 emission, however, is found 
to have an emission region of size 0s Z 10". Although the latter 
violates the original assumption that 0, <K 0beam, it does indicate 
that CO 7-6 stems from a much more extended portion of the 
envelope than HCN 9-8. Implications of the sizes of emission 
regions derived here are discussed in Sect. 6.4 



4. Envelope spectral line energy distribution model 

4. 1 . Model setup 

Radiative transfer calculations with the Monte Carlo ratran 
code (Hogerheijde & Van der Tak 2000) are employed to pro- 



Line transition 


measured T^b (K) 


Tnib ratio 


inferred 0s 




JCMT Herschel 




(") 


CO 7-6 


88±44 22±2 


4±2 


>10 


HCN 9-8 


9±4.5 0.32±0.03 


28+;^ 
>10 


«1 


HCN (V2= 1)9-8 


4+2 <0.2 


<3 



Notes. JCMT intensities are taken from|Boonman et aLH2001f, Her- 
schel intensities are from this work. The HCN (v2 = 1) 9-8 upper limit 
for Herschel is derived from the ~0.1 K l-cr noise level near 797.3 GHz 
in our spectrum. The inferred size of the emission region, 0„ is calcu- 
lated using Eq. |TJ. 



the LAMDA database (Schoier et al.| |2005 
vided by Yang et al 



which uses data pro 
(2010) for CO, IDumouchel et al.' (2UW. 



>0T0| 
et aL\ 



3.2 



vide a model comparison to the observed SLEDs from Sect. 
We rely on spectroscopic parameters and collisional rates from 



for HC N and HNC, .Flower, (1999) for HCO+ and [Turner 
([1992) for CS. 

The physical model we use is derived from a fit to the ra- 
dial intensity profile of the sectors of the 450 and 850 jum JCMT 
SCUBA maps that appear spherical and to the /I > 50 /urn sec- 
tion of the continuum spectral energy distribution (Luis Chavar- 
ria, private communication), with a strategy similar to that of [Van[ 
[der Tak et al.[ ( [1999[[2000I ). Our current model takes into account 
the new distance of 3.3 kpc and a corresponding input luminos- 
ity of 2.1x10'' Lq. Briefly, the model consists of a radial density 
gradient for molecular hydrogen described by an r ' power law 
(4x10^ down to 7x10- cm"-'), a total gas mass of ~500 Mq, and 
a corresponding temperature profile which varies monotonically 
from ~ 1400 K in the innermost shell to 27 K in the outer shell at 
70 kAU. The radial distance at which the signal in the SCUBA 
450 fim map falls below 3cr determines the value of the outer ra- 
dius of the model envelope. The total luminosity, as determined 
by integrating over the best fit model continuum spectral energy 
distribution, is 8x10"* Lq; the difference with the observed lumi- 
nosity is due to the absence of a cavity and a disk in our one- 
dimensional model, which does not include near-infrared emis- 
sion. jJimenez-Serra et al. (20121 have found chemical segrega- 
tion at very small scales (<!", 3 kAU), but this is not expected 
to impact our models at the scales probed by our single dish ob- 
servations (>15", 50 kAU). Therefore, molecular abundances in 
our models are kept constant throughout the envelope and are 
scaled to match the observed SLEDs. 

After the Monte Carlo method solves for the molecular level 
populations, ray tracing determines the emergent line intensity, 
which is convolved with a Gaussian telescope beam with the ap- 
propriate HPBW at the relevant frequency. Resulting spectral 
maps are collapsed along the spectral direction to obtain the in- 
tegrated line intensity value, which is compared to observations 
in Fig. [4] Because the outflow component is not represented in 
the model, we compare predicted integrated line intensities from 
the models described above only with the envelope component 
of the observed spectral lines (Sect. 3.2 1. 

Despite the obvious increase in linear size scales when com- 
paring an older model for AFGL2591 assuming a distance of 
1 kpc (e.g., [Van der Tak et al.||2000J [Van der Wiel et al.|[20TT) 
with the one employed in this work, the differences in emergent 
line intensities are relatively minor. This is due to the similar 
temperature and density structure, which lead to comparable ex- 
citation conditions. A noteworthy element of the new model is 
the temperature of >1000 K in the central shells of the envelope, 
compared to 372 K for the old model. The maximum tempera- 
ture in the model is essentially set by the arbitrary choices of the 
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thickness of the spherical model shells and the size of a central 
sphere which is vacant of gas. By running test models in which 
the very hottest inner shells are removed, we determine that the 
hot gas in the central envelope has no impact whatsoever on the 
emergent line intensities of rotational levels up to £'iip=700K. 
Essentially, none of the Une emission emanating from the hot 
gas escapes the envelope to be detected by external observers. 

4.2. Comparison of model results with observations 

With a CO abundance of 5x10"^ rel ative to H?, one quarter o f 
the 2x10"** applicable for warm gas (' Van der Tak et al.||l999| ), 
the modeled SLED of the main form of CO Ues a factor of 2 to 4 
above the observed points up to 7=9-8 (off the scale in Fig. [4^). 
Beyond that, the model SLED drops off sharply and underpro- 
duces observed line intensities from 7=13-12 onward by factors 
of 2 to 3. Although with smaller discrepancies, the models for 
the isotopic variants '^CO, C'**0 and C'^O (with abundances of 
5x10"^, 6x10"*^ and 3x10"^, respectively) show a similar trend 
in Fig. |4|3,c,d: low-J emission is over-predicted while high-7 
emission is under-predicted. 

The HCO^ model, at an abundance of 1x10"^, seems suc- 
cessful: it matches the envelope component of the observed lines 
7=6-5 and higher to within a factor 2 (Fig. |4^). However, the 
same model overproduces HCO^ 4-3 by a factor 5. The mod- 
els for the other species studied here, CS, HCN and HNC (abun- 
dances of 4, 6 and 1x1 0"*^ , respectively), again overproduce low- 
7 emission by factors ~5 and are marginally consistent with or 
fall below the observed points at higher 7 (Fig.|4j',g,h). 

The general picture that emerges is that our smooth, spher- 
ically symmetric models predict too much emission from cold 
gas (~50K) or too little emission from warm gas (>100K). This 
effect persists even when we consider the possibility that the en- 
velope contribution to the emission is overestimated by a fac- 
tor of 2 for lines that are fitted by a single Gaussian (gray error 
bars in Fig. |4]i. For all species studied here, choosing a differ- 
ent molecular abundance has an effect on the vertical offset of 
the model SLED, especially at more optically thin transitions at 
high energies, but it does not sufficiently reduce the slope be- 
tween £'up~100 and 200 K to match the observations. Our model 
serves to show that a spherical, homogeneous geometry offers a 
poor match to the true excitation balance. Therefore, the abso- 
lute molecular abundances used in this section are not meaning- 
ful by themselves. For example, the CO abundance in our model 
is lower than generally assumed, but stiU severely overproduces 
low-7 line intensities and optical depths. 

A model might be considered where CO is depleted from 
the gas phase, to suppress specifically the low-7 CO emission. 
However, the temperature profile we use for the envelope does 
not reach values below 25 K. Hence, the model gas is too warm 
for CO freeze-out to occur, consistent with the absence of CO 
ice features in mid-infrared spectra ( [Van der Tak et al.|1999 1. 

Our model is based on a density profile with a power law 
index of a = -1.0, as outlined in Sect. |4.1| and also used in ear- 
lier work on this source (see Sect.[T]l. We have explored ratran 
models with steeper profiles (a = - 1 .5), but the resulting SLEDs 
are very similar in shape to the ones with a - -l.O. A steeper 
power law therefore does not mitigate the mismatch of the shape 
of the linear rotor SLEDs with the observations. 

Considering the spectral dimension, the modeled and beam- 
convolved line profiles of, e.g., CO and HCO^ up to 7-6 and 
HCN up to 7=10-9 appear broadened or even self-absorbed, 
indicating optically thick emission. For example, the optical 
depths at line center determined during the ray tracing, for the 



main isotopic form of CO, are between 90 and 10 for 3-2 up 
to 11-10, while only the 16-15 is marginally optically thin at 
0.8. The modeled lines for the isotopologues of CO are all op- 
tically thin, as are all CS transitions except 7=7-6. However, 
HCO+ is optically thick until 10-9, HCN until 9-8 and HNC 
until 7-6. Indications for such high optical depths are not seen 
in the observed line profiles (Figs. [T and|2]i. This discrepancy 
between observed and modeled line shapes confirms earlier sus- 
picions ( |Van der Tak et a l. 1999; Van der Wiel et al. 12011 ) that 
line optical depths are too high in a spherical model. In Sect. |6.1| 
we will suggest possible approaches to address the optical depth 
with more sophisticated models. 



5. Physical conditions per component 

5.1. Foreground clouds 

Our CO 5-4 spectrum exhibits an absorption feature near 
Okms-i (Fig.|7]l, which is also seen in JCMT data of CO 3-2 and 
is coincident in velocity space with a known foreground cloud 
( |Van der Tak et al.|I999| l. Conversely, an emission peak is visible 
near Okms"' in our IR AM spectra of '^CO I -O and C"^0 1-0 
(Eup = 5.3 K, Fig. |7]l. [ Mitchell et aL| ( p^92| show spectra of 
CO 2-1 toward various positions m the envelope showing the 
same absorption feature near Visr = Okms"'. It is not seen in 
higher-7 lines of CO, indicating that the absorbing foreground 
cloud must be relatively cold. It must also be of relatively low 
column density, given that the absorption feature is not present 
in i ^CO and C'^O 5-4. 

Minh & Yang (20081 observe an extended Okms ' compo- 
nent in CO and '^CO 1-0 emission, peaking 8' to the northeast 
and northwest of AFGL2591. Our observations show that this 



gas cloud also extends to the position of AFGL2591 itself, and 
that it must be in the foreground, since it absorbs CO 3-2 and 
5-4 emission from the warm protostellar envelope. 

We construct a simple 'slab' component to account for the 
foreground CO absorption, using the interface to radex ( |Van| 
der Tak et al. 12007 ) available in cassis. First, the shape of 



the CO 5-4 emission line is mimicked by three 'LTE' com- 
ponents with excitation temperatures of 100 K, Vcentr of -4.5, 
-7 .3 and -I- 1 .6 km s- 1 , FWHM values of 5 . 1 , 1 5 .7 and 2.5 km s" ' 
and CO column densities of 6, 9 and 0.3x10'^ cm"^. Note that 
these emission components are merely invoked to create repre- 
sentative background emission and are not meant to represent 
real physical conditions in the protostellar envelope. On top 
of the mock background emission, a radex slab is added with 
A^co = 3x10'^ cm"^ and a kinetic gas temperature (r^n) of 10 K. 
This choice for temperature is motivated by the fact that it should 
be lower than the minimum gas temperatures of ~25 K reached 
at the outer edges of the protostellar envelope of AFGL2591, 
while at the same time it should be high enough to result in de- 
tectable emission in the 1-0 transition of '^CO and C'**0, with 
£'up=5.3 K. 

Because the foreground component is likely to be part of a 
larger scale cloud ( Minh & Yang|2008) , the angular size of the 
slab is set at 50", larger than any of the beam sizes of our ob- 
servations. Finally, the absorption component is given an intrin- 
sic FWHM of 0.8 km s"' and V^i of +0.3 km s"' to match the 
width and position of the observed absorption feature. After the 
frequency-dependent optical depth calculation and the radiative 
transfer, the result is convolved with the relevant telescope beam 
to make a direct comparison to observed line profiles in Fig. [7] 

With the parameters for the cold material chosen above to 
match CO 5-4, we also find self-absorption in CO 3-2 (again. 
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Fig. 7. Zoom of the observed km s"' feature in various CO lines 
(histograms). The results from the radex slab are shown as solid lines. 
To improve visibility, the brightnesses for C'*0 (both observations and 
models) are multiplied by 10. Vertical offsets are applied as follows: 
+5 K for CO 5-4, -5 K for C'**0 2-1, -10 K for C'**0 1-0 and -15 K 
for '^^CO 1-0. The gray dashed line marks the approximate centroid 
velocity of the absorption feature at +0.2kms"'. The model represen- 
tations only include a component for the main envelope emission for 
5-4 and 3-2, where a background is needed to absorb against. The 
2-1 and 1-0 model lines only account for the foreground emission at 
~Okms"', hence the discrepancies at V^r < -1 kms"' where the main 
envelope of AFGL2591 dominates. 



against a mock background emission line) that is qualitatively 
consistent with the JCMT observations. Moreover, due to the 
10 K gas temperature, the same component naturally explains 
the emission bumps observed in '^CO 1-0, C'^O 1-0 and even 
in C'**0 2-1 (Fig.pb. The column densities needed for '^'CO and 
C'^O are 3x10'^ and 2x10'** cm"^, respectively. 

We conclude that a cold slab of foreground gas explains the 
0km s"' feature. Assuming an abundance of CO/H2 of at most 



10 , we derive a column density A^h, ^ 3 x 10" cm" . The 
visual extinction (Ay) for this column density would be ~3 mag- 
nitudes, following |Bohlin et al.| ( |1978| ), purely for the molecular 
part. If we add to this 1 magnitude each (front and back) for the 
atomic surfaces of the cloud where CO is photodissociated, we 
find that the total Ay > 5. The derived conditions in the km s"' 
component ('foreground A') are summarized in Table [4] 

Whereas the 1-0 emission features are insensitive to changes 
in Tkin between 5 and 20 K, the depth and shape of the 
CO 5^ absorption features constrains to [9,15] K. The 
value adopted for A^co is motivated by the depth of the self- 
absorption feature, which starts to deviate significantly when 
A^co is changed by a factor > 2. The model result is insensitive 
to the number density of the collision partner in radex as long 



In fact, if any feature is present at that velocity, it is in emis- 
sion (Fig. [T]i. The slab model with r^in = lOK and A^co = 
3 X 10'^ cm"^, however, predicts a detectable absorption feature 
for CO 6-5. We are unable to explain the CO 6-5 without com- 
promising the match to the lower-7 0km s"' feature. 

Although not included in the considerations in this paper, 
several additional probes of the Okms ' foreground cloud are 
available. Firstly, besides the CO signatures shown in Fig.|7] the 
^Pi-^o line of neutral carbon (at 492.161 GHz in our CHESS 
spectral scan) also shows a weak emission feature near km s"' . 
Secondly, an upper limit on the density might be derived from 
the fact that the Okms ' component is not seen in lines of, 
e.g., HCO+ and HCN (Fig. |2|. Thirdly, signatures of CH, 
CH^, H^ O^ and other hydride sp ecies do show absorption near 



Okms ' d Bruderer et al.||201 0b), as does the HF signature in 
the CHESS spectrum detected by Emprechtinger et al. (2012). 
A more extensive exploration incorporating all observed tracers 
of the Okms ' component toward AFGL2591 is therefore war- 
ranted. 

In addition to the km s ' component described above, there 
is another absorption component, at H-13kms visible in Her- 
schel observations of HF (Emprechtinger et al. 2012, but not dis- 
cussed there) as well as H2O ( Choi et al.,in preparation) . It is not 
seen, however, in any lines of CO species presented in this work, 
neither those observed with HIFI (Fig. [TJ nor those observed at 
IRAM 30m and JCMT (the velocity range covered in the obser- 
vations is larger than what is shown in Fig.|7]i. The diatomic hy- 
drides observed in the Herschel WISH program ( [Bruderer et al. 
|20I0b| l do not show clear evidence of emission or absorption at 
-1-13 km s"' either. 

Based on our '^CO 1-0 spectrum, the -i-13kms ' feature 
would need to have a peak intensity of <0. 1 K to ensure a 
non-detection. With that upper limit, a brief exploration of 
RADEX models tells us that the '^CO column density is at most 
10'"^ cm"^. Using the canonical scaling factors to convert to CO 
(60) and noting that the CO abundance relative to H2 can be as 
low as 10""^ in diffuse gas, we derive that Ny^-, < 6x10^' cm"^ 
in the H-13kms ' component. In addition, the CHESS spectrum 
of HF 1-0 presented by Empr echtinger et aL] p012) provides 
a low er limit for the coluri m density of the -1-13 km s ' cloud. 

do not discuss the -1-13 km s"' com- 



Emprechtinger et al. 



40^ 



cm , i.e., when the low-7 CO excitation (e.g., iChoi et al. 



as «(H2) is above • 
balance is in LTE. 

Given the inferred nature of the foreground cloud, it is im- 
portant to consider that H2 densities in translucent clouds are 
generally in the order of 10^^ rather than > 10"* cm"^. With a lower 
H2 density, the 7=1-0 emission features in '^CO and C'^O can 
be explained with temperatures as high as 17 K. Hence, the den- 
sity as well as the upper limit to the kinetic temperature of the 
foreground cloud are poorly constrained. 

In addition, a discrepancy is found with respect to the ob- 
served CO 6-5 line profile that shows no absorption at km s ' . 



While 

ponent, an absorption feature does appear at that velocity with 
a depth of ~70% relative to the saturated Okms ' component. 
For the latter, they derive A^hf > 4x10'"' cm"^. Scaling by the 
relative depth of the absorption features and assuming that all 
fluorine is in the form of HE (HF/H2=3.6xlO"^), we find that 
A^H. ^ 8x10^" cm"^. Column density limits are summarized in 
Table |4]( 'foreground B'). 



5.2. Outflow gas 

Previous studies indicate that CO species, HCO^, and neutral 
atomic carbon are present in the outflow in significant amounts 
[I994t pasegawa & Mitchell. 1 995, |Doty et al.| 



2002). Our results confirm that line profiles of CO species and 
of HCO^ are well represented by two separate kinematic com- 
ponents (Sect. [3^. 

Using radiative transfer calculations with the radex code, we 
estimate kinetic temperatures, number densities and column den- 
sities of CO and HCO^ in the outflow component. This approach 
requires an assumption of an isothermal homogeneous material, 
in which the radiative transfer is solved using collisional rates 
from the LAMDA database (|Schoier et al.||2005||. Following 



|Van der Tak et al.] ( |2007[ ), we construct grids of radex models 
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Table 4. Physical conditions separated by component. 



Component 


Vlsr 


FWHM 






Nco 




X(HCO+)^ 




(kms-1) 


(kms-') 


(K) 


(cm-3) 


(cm^^) 


(cm^^) 




Foreground A 


+0.2 


0.8 


lO^i' 




>3xlO" 


>3xlO^' 




Foreground B 


+ 13 








<6xlOi5 


-8x1020-6x10^' 




Outflow 


-7 


-6-10 


60-200 


-10^-10^ 


1.5x10"*'' 


1.5x10^2 


-10-9-10-^ 


Envelope"^ 


-5.7 


-3-4 


23-372 


-lOMxlO^ 


-10'** 


-10^2 


4xl0-'"-lxl0-'* 



Notes. A '-' symbol indicates that the corresponding quantity is unknown for that component. 

Molecular abundance of HCO^ relative to H2. 
C'* Scaled from A^c'^o = 3 x 10'^ cm^^. 

Envelope number density and kinetic temperature ranges are taken from the parameterization of the spherically symmetric model 
(Sect.|4li. 




CISq 5-4/7-6 
HCO+ 4-3/7-6 



3 4 5 6 7 8 

log[«H2] (cm"-'') 

Fig. 8. Lines of constant intensity ratio of C"*0 5-4/7-6 (solid lines) 
and HCO* 4—3/7-6 (dashed lines), resulting from a radex grid calcula- 
tion. Curves are labeled by their line ratio value. The gray areas mark 
the l-o" uncertainty ranges of the observed line intensity ratios in the 
outflow component, including the fit error (see Table [2j and the 10% 
calibration uncertainty (see Sect.jZTJ. The darker gray area denotes the 
region where the observed line intensity ratios of both tracers overlap. 



and calculate line intensity ratios for values of between 10^ 
and 10** cm"^ and r^in between 10 and 400 K. In the input for 
the grid calculation, the column densities are intentionally low 
in order to keep all modeled lines optically thin. In this limit, 
the intensity of each line is directly proportional to the column 
density and line ratios are therefore insensitive to the absolute 
column density. In Fig. [8] we show the resulting line ratios of 
HCO+ 4-3/7-6 and C"*0 5-4/7-6 as a function of riu, and Tki,,. 
When the observed line intensity ratio is indicated in the dia- 
gram, it is possible to constrain a region in parameter space that 
fits the observations. 

First, the mere presence of high-/ HCO^ emission in the out- 
flow component is an indication that number densities must be 
relatively high, given that critical densities are in the order of 
10^ cm"^ for the 6-5, 7-6 and 8-7 transitions. Second, from 
the C'^O line ratio in Fig. [s] it is evident that densities below 
lO'^cm"^ require temperatures that are unphysically high. The 
observed HCO^ and C'**0 line ratios intersect in a region in 
(nHT,7^kin)-space near lO'^'^-lO^ cm"^ and -60-200 K. Inside this 
region, higher densities imply lower temperatures, and vice versa 
(Fig. [sj. The measured line ratio for '^CO 5^/7-6 gives the 
same constraints as C'^O, confirming that signal-to-noise and 



optical depth do not affect our result. Finally, a line ratio diagram 
for 1-0/2-1 suggests the presence of an additional compo- 
nent of colder gas at Tkin < 50 K. However, this result should 
be considered with caution, because (i) the kinematic separa- 
tion of envelope and outflow contributions is not unambiguous 
(Sect. 3.2 1, and (ii) the beam filling factor of the outflow compo- 
nent could be different for the 2-1 and 1-0 observations (beam 
HPBW 11" and 22", respectively). 

With «H,=[10^10^]cm-2 and Tyn = [60,200] K, a column 
density in the order of 10'^ to lO'^'cm ^ for HCO^ reproduces 
the line strengths observed for the outflow. Similarly, a C'^O 
column density of 3x10'^ cm"^ (+30%) is required to explain 
the HIFI C"*0 line strengths in the outflow component. With a 
CO/C'^O ratio of 500 and a CO abundance of 10 relative to 
H2, we derive A^h, = 1.5x102^ cm ^ and a molecular abundance 
in the order of 10"' to 10"^ for HCO^ in the outflow gas. This 
HCO^ abundance falls in the same range as pr eviously adopted 
values for th e envelope, ranging from 4x10"'° (Carr et al. 1995 1 
to -1x10-^ dVan der Wiel et al.|201l| . 

In conclusion. Table |4] contrasts the physical conditions de- 
rived here for the protostellar outflow component (Sect. |5.2] i with 
those found for the CO foreground cloud (Sect. |5.1| l and in the 
envelope. Our current analysis indicates that the outflow gas is 
not significantly different from that in the envelope, considering 
gas density, gas temperature, as well as the chemical balance of 
CO and HCO+. 



6. Discussion and conclusions 

6.1. Envelope structure 

The spherically symmetric radiative transfer models are unable 
to simultaneously explain the low- and high-7 transitions of the 
observed integrated line intensities (Sect.|4] Fig.|4]i: the modeled 
SLEDs of all species studied here decline too steeply between 
-100 and 200 K. This indicates that the passively heated, spher- 
ical envelope structure does not populate the high J levels suffi- 
ciently, or that emission from those levels is too optically thick. 
Line optical depths are known to be higher than observed for 
AFGL2591 in ground-based low excitation lines ( |Van der Wiel| 
|et al.|20l T) and in Sect. [4] we have noted that optical depths for 
the modeled higher J lines are lower, but still >1 in many cases. 

To boost line intensities at E^p above -150K without in- 
creasing the lower level intensities too much, it would be ap- 
propriate to apply a model geometry that includes outflow cav- 
ities. This would affect not only the density structure to create 
lower opacity pathways (if a suitable inclination angle is cho- 
sen), but also the chemical and thermal balance through direct 
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illumination of envelope material by UV radiation from the cen- 
tral massive protostar. Indeed, we demonstrate in Sect. [4] that 
the hot gas (>1000K) in the center of the model envelope does 
not contribute to emergent line intensity if it is embedded in an 
isotropic envelope. 

In this hght, our AFGL2591 observations should be com- 
pared with a model scheme such as that by |Bruderer et al.| 
(2009a 2 010a|[2012| l, with a physical geometry appropriate for 
the massive protostellar envelope of AFGL2591. Their model 
allows for an outflow cavity to be present in an axisymmetric 
(or even three-dimensional) density distribution, from which the 
dust temperature profile is derived through continuum radiative 
transfer. The UV field is calculated and used as input to itera- 
tively determine the position-dependent chemical balance, ex- 
citation balance, gas temperature and resulting line emission. 
The fully self-consistent chemical treatment includes freeze-out, 
photodissociation. X-ray ionization and other processes. Prelim- 
inary exploration of similar models already shows that, as the 
cavity allows a larger volume of lower density gas to be irradi- 
ated and heated directly, SLED shapes become significantly flat- 
ter compared with spherically symmetric geometries (Van der 
Wiel|201 1 , Chapter 4; Steven Doty 2011, private communica- 



6.2. Foreground clouds 



tion). Bruderer et al.| ( [20T2| also show that in exposed layers of 
the gas envelope, the gas and dust temperature become decou- 
pled, rendering inadequate our model for which the gas temper- 
ature is derived from dust continuum measurements (Sect.|4|. 

The HCO^ molecule in particular could provide an interest- 
ing probe of outflow cavity geometry, since, contrary to CO 
and HCN, the HCO^ abundance is expected to be enhanced 
in the cavity walls (Simon Bruderer 2011, private communica- 
tion). Additional constraints could come from the incorporation 
of CO^ observations. Although no signatures from this species 
are detected in the spectral scans with HerschelfHlFl ( Paper II| l 
and with JCMT/HARP-B ( |Van der WielHIOTT] Sect. 3 .A ), two 
submillimeter CO^ lines are detected in deeper spectra by'Stau-J 
|ber et al.] ( i200 7 ). From the sample of low-mass and high-mass 
protostars in their study, they conclude that X-rays play a role 
in the dissociation and excitation of molecules in low-mass pro- 
tostellar envelopes, but that CO^ observations can be explained 
by far-UV radiation in high-mass objects. Even in dissociation 
regions, CO^ emission could be weak or absent, since this ion is 
easily destroyed by collisions ( [Stauber & Bruderer|2009[ ). 

Additionally, line opacity properties of the modeled proto- 
stellar envelope can be adjusted by distributing the same amount 
of gaseous material inhomogeneously, as previously suggested 
by |Van der Wiel et al.| ( |201 1") . This could be used to further alle- 
viate the optical depth effects (Sect.|4]l, if the introduction of cav- 
ities in a two-dimensional geometry does not prove sufficiently 
effective to match the observations. 

Studies of low-mass star-forming envelopes have shown that 
a passively heated envelope is the dominant component for lines 
with < 200 K, but that the -200^000 K section of the CO 
SLED is a discriminating probe of the contribution of UV ir- 
radiation and shocks (Van Kempen et al.||2010[ |Visser et alj 
|2012). While the highest / CO line presented in this work is 
16-15, the highest attainable to HIFI, Herschel/PACS has ob- 
served AFGL2591 at frequencies covering lines up to /up ~ 30, 
to be analyzed in the context of the WISH key program ( |Van| 
Dishoeck et al.|20TT] l. We expect the extension of the CO lad- 



We identify a Okms ' foreground cloud in Sect. 5.1 for which 
we derive a kinetic gas temperature between 9 and ~17 K, based 
on the combination of absorption and emission signatures in the 
J-l-0 up to 5^ transitions of CO and isotopologues. The ma- 
terial in the foreground cloud discussed here is redshifted with 
respect to AFGL2591. Therefore, the foreground could be feed- 
ing additional gas to the protostellar system, but evidence would 
have to be found to support physical proximity of the Okms ' 
cloud to AFGL2591. 

For a second foreground cloud, at Visi— -h13 kms"', we use 
our non-detections in CO and the absorption in HE to derive lim- 
its for the column density: 8x10^" cm"^ < A^h, ^ 6x10^' cm"^ 



(Sect. 5.1 



The CO abundance relative to H2 must be < 10"^ 
for the derived limits to be consistent with each other (for 
CO/H2=10 the upper Hmit would be 6xl0''cm"^). In such 
diffuse conditions HE is more abundant than CO, as predicted 
by N eufeld et al.|(|2005) l and observed, for example, by |Sonnei>] 
itrucker et al. (2010 1 on the line of sight to W51. It is surpris- 
ing, however, to detect a relatively strong absorption signature 
of H2O without the presence of CO. Water abundance values 
will be derived by |Choi et al.| ( |in preparation] ). 



6.3. Outflow gas 

The kinetic temperature of 60-200 K and number density of 
^lO^-lO*" cm"^ derived for the outffow gas are comparable to 
those in the passive envelope (Sect. |5.2[ ). The component that 
we identify a few kms"' blueward of the systemic V]sr is kine- 
matically distinct from the envelope (Table [2]), but it apparently 
harbors gas that is similar in density and temperature. Our re- 
sults are consistent with conditions found in other protostellar 



outffows (e.g., Giannini et al.||2001| 


Moro-Martm et al. 


20011 


Lefloch et al.|2010||Yildiz et al.|2010| 


2012t|Bjerkehetal. 


2011 


Gomez-Ruiz et al.||2012|l. Specifically, CHESS spectra toward 



10^cm~^ and temperatures around 200 K (jCodella et al.||2012 



[Benedettini et al.|2012[|Lefloch et al .120101120121 1. 

We speculate that the similar physical and chemical condi- 
tions (Table]?]) between outflow and envelope could point to out- 
flow gas that has been expelled from the envelope relatively re- 
cently. The blue component observed in this work may in fact 
trace swept-up envelope material (<10kms"' from the systemic 
velocity) rather than fast outflowing gas that is farther separated 
from the parent envelope, see, e.g.. Van der Tak et al. <\1999) for 
AFGL2591, or jQiu et al.l ( ]20TT] ) and jGomez-Ruiz et al.l ( 120T2] i 
for other protostars. Our outflow conditions are compatible with 
the temperature and density derived for the warmest of the three 
outflow components in L1157-B1, as found by Le floch et al.] 
( |2012 ! ). These authors argue that lower kinetic temperatures indi- 
cate older outflows, analogous to our conclusion that the warm, 
dense outflow material from AFGL2591 is young. We should 
also caution that the double Gaussian fits that are applied to 



der to higher energies to shed more light on excitation drivers in 
massive protostellar envelopes such as AFGL2591. 



distinguish the components (Sect. 3.2 1 are likely to suffer from 
cross-contamination. Both components could partly trace the 
same gas and, as a consequence, the derived physical conditions 
for the envelope and outffow may be artificially equalized. 

6.4. Compact HCN emission 

In Sect. |3.4[ our observations of HCN 9-8 and CO 7-6 are com- 
pared with those from a larger telescope, in order to impose lim- 
its on the size of the emitting regions. We find sizes of <1" for 
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HCN 9-8 and >10" for CO 7-6. |Benz et al.| ( |2007] l find source 
sizes of <l" in an interferometric HCN 4—3 map, but the ex- 
tended emission in their observation is likely filtered out. Con- 
versely, single dish observations by | Van der Wiel et al.| ( |20lT) 
reveal larger scale (>30") emission in the 4-3 transition of HCN, 
but this line traces lower temperatures and lower densities than 
HCN 9-8. Considering our envelope model from Sect. |4] only 
in the very inner shells (< 30 AU ~ 0.01") does the number den- 
sity come within an order of magnitude of the critical density 
of HCN 9-8, 3x10*^ cm"^. It is therefore conceivable, and fully 
consistent with our observations, that the HCN 7=9-8 emission 
arises in a region smaller than 1" in angular size. The critical 
density for CO 7-6 is a factor -50 lower than that of HCN 9-8; 
it is therefore expected that CO 7-6 emerges from a more ex- 
tended region of the envelope. 



7. Summary and outlook 

In this first Herschel CHESS paper on the massive star-forming 
envelope AFGL2591, we present HIFI spectral line measure- 
ments of a wide range of transitions (Eu^jk-^Q up to 300-700 K) 
of CO, '^CO, C'^O, C'^O, HCO+, CS, HCN and HNC, sup- 
ported by ground-based measurement of lower excitation lev- 
els. Emission lines are separated into a narrow 'envelope' and a 
broader 'outflow' component, if the S IN is sufficient (Sect. [T2) i. 
The integrated line intensities of the envelope component, pre- 
sented in the form of spectral line energy distributions (SLEDs), 
are compared with radiative transfer results based on a spheri- 
cally symmetric envelope model (Sect. |4]). 

The main conclusions of this work are the following. 

1 . The model SLEDs for all species fall ofi^ too steeply between 
100 and 200 K to explain the flat shape of the observed en- 
velope SLEDs (Sect. |4]i. We expect that a model including 
outflow cavities and direct UV illumination will provide a 
better match (Sect. [6T| . 

2. The outflow gas is found to be as dense and warm as that 
in the quiescent envelope (Sect. |5.2| l. We conclude that the 
outflow gas traced here may have been expelled only recently 
(Sect.[0]l. 

3. Trends of line centroid and line width with frequency, ex- 
citation level or telescope beam size are identified in the 
envelope component. In absence of a single, consistent, 
physical explanation, we ascribe these trends to ambiguities 
in the separation of the envelope and outflow components 
(Sect.[33]l. 

4. By comparing line intensities measured with Herschel 
(3.5 m) and JCMT (15 m), we conclude that HCN 9-8 
emission must originate from a region <1" in angular size 
(Sects. [141 and|64ll. 

5. A known foreground cloud at Visr^Okms"' is detected in 
several lines of CO and isotopologues. Limits to its temper- 
ature and column density are derived in Sect. 5.1 

6. A second foreground cloud, detected in HE and H2O 
at -1-13 km s"', remains undetected in our observations 



(Sect . \5.\) . Its nature and column density are discussed in 
Sect.l62r 



Future papers in this series based on the CHESS observations 
of AFGL2591 will be ordered by chemical families of species, 
e.g., N-bearing species (Paper II), S-bearing species. |Paper II 
will also include details of all other molecular line signal de- 
tected in the 490-1240 GHz spectral scan. In addition, future 



work should be devoted to the application of a more sophisti- 
cated model geometry for AFGL2591 to explain the highly ex- 
cited rotational lines accessible with Herschel. Such a model 
should include at the least an outflow cavity geometry and chem- 
ical effects of direct UV illumi nation (n ow ava ilable for hydrides 
and C/C+/CO, [Bruderer et al. "' 2010a b , | 20T2| , but possibly also 
shock physics and inhomogeneous envelope structure. 
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